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We propose a simple structure for stationary non-Markovian quantum chains in the framework of 
collisional dynamics of open quantum systems. To this end, we modify the microscopic Markovian 
system-reservoir model, consider multiple collisions with each of the molecules with an overlap 
between the collisional time intervals. We show how the equivalent Markovian quantum chain 
can be constructed with the addition of satellite quantum memory to the system. We distinguish 
quantum from classical non-Markovianity. Moreover, we define the counts of non-Markovianity 
by the required number of satellite qubits and bits, respectively. As the particular measure of 
quantum non-Markovianity the discord of the satellite w.r.t. the system is suggested. Simplest 
qubit realizations are discussed, and the significance for real system-environment dynamics is also 
pointed out. 



I. INTRODUCTION 



In recent years we have seen increasing interest in abstract modeling of time evolution of open quantum systems 
coupled to their environment via different quantum channels. While decrypting fundamental laws is their ultimate 
goal, inspirations of such models and possible future applications of their results range from solid state physics to 
quantum biology and quantum information technology [IH1]. 

Dynamics modeled microscopically as a series of discrete interactions or "collisions" between a central system and 
environmental (reservoir) molecules leads to abstract collisional system-reservoir models. For a memoryless Markovian 
time evolution the reservoir has the idealized ability of forgetting new information faster than the rate of collisions; 
i.e., after each collision the reservoir would totally relax to its pre-collisional state before the next collision occurs. 
For the (central) system, this results in what we call Markovian quantum chain. In each step along the chain, the 
system undergoes decoherence, leakage of information out to the reservoir to be forgotten there. Mathematically, 
the Markov quantum chain's irreversible dynamics is obtained from reduction of the unitary time evolution of the 
system-environment compound by means of partial trace over degrees of freedom of the environment. In the special 
case of stationarity, the mathematical tool of semigroups of completely positive (CP) trace preserving dynamical maps 
can be used [2H6]. 

The collisional model of Markovian chains has been studied in numerous works [7J [5] . It will be particularly useful 
if we want to monitor the system [51-11 1| . Although we are never allowed to directly measure the system state, we can 
measure each reservoir molecule after its collision with the system. By their regular selective measurements we can 
monitor the system state indirectly. The corresponding resolved evolution of the quantum chain is called selective 
Markovian quantum chain. 

A more realistic description requires an account of the reservoir memory and the induced non-Markovianity of the 
dynamics of the central system (the chain). For analyzing new dynamical properties and quantifying the difference 
from the simple Markovian processes, several non-Markovianity measures have been proposed so far. These make use 
of different aspects of non-Markovian (NM) evolutions, e.g., the non-divisibility of the underlying quantum dynamical 
map |12j . the increasing trace distance (i.e.: distinguishability) of two initial quantum states, accompanied by the 
back-flow of information from the reservoir to the system [13], or the discord between system and reservoir [14j . 
Non-Markovianity is a field of active debates (cf., e.g., [T5]) so we emphasize the importance of abstract models to 
capture fundamental aspects of memory keeping processes. 

In this paper we study the structural features of NM open quantum systems by constructing discretized NM 
processes, i.e., quantum chains. Starting from the generic structure of a Markovian quantum chain we are going to 
impose a certain NM structure on it. The proposed abstract NM structure corresponds to an open quantum system 
in a reservoir of non-interacting molecules which only collide with the central system through unitary collisions. The 
initial, for the moment Markovian, quantum chain can be seen in Fig. 1(a) To engineer the NM element, we allow 
each molecule to collide with the system more than a single time, say twice, while the molecules' collisional intervals 
overlap, by assumption, with the collisional periods of both previous and consecutive molecules, cf. Fig. 1(b) Instead 



of separating what is forgotten and what is not by different timescales, our abstract model has a built-in exact 
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FIG. 1. Quantum chains: different system-reservoir collision models (g: 
molecule unitary collision operator), (a) Standard Markovian ch 
once, (b) The proposed NM ch ain: the system interacts with each independent molecule twice, with overlap between collisional 

(c) NM chain in [18] : the system interacts with each molecule only once, but the molecules are 
the system interacts with the same molecule multiple times. 



intervals of 'nearest' molecules 
initially entangled 



(d) NM chain in [19] 



memory time. Most importantly, we show how our NM chains become Markovian at the price of attaching satellite 
memory subsystems. Due to our proposed specific NM structure, the identification of the satellite memory within the 
environment is straightforward compared especially to the efforts requested in oscillatory reservoir models (cf. [161 117] 
and Refs. therein). Then a natural count of non-Markovianity follows: let it be the size (in qubits) of the requested 
satellite memory. We can refine this count into informatic measures of non-Markovianity, as we show later. 



1(d) Ref. [IS] introduces NM mechanism by starting all molecules from an a priori entangled state, cf. Fig 
work focuses on the molecular realization of a single NM step instead of general NM chains. In [ 
consists of a single molecule and this molecule collides with the system consecutively, many times, cf. 
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important advantage of our NM structure over these two is that it can invariably host the monitoring of the system, 
which was so instrumental for quantum Markov chains but becomes complicated in structure Fig. 1(c) and impossible 
in Fig. |l(d)| 

In Sec. [II] the exact notion of a non-selective and a selective quantum Markovian chain is defined. The definition is 
extended for quantum NM chains in Sec. |III[ Then Sec. |IV| contains our new results demonstrated on the toy models 
of one-qubit NM chains. Based on these results we discuss in Sec.|V]a new measure of quantum non-Markovianity as 
well as a distinction between classical and truly quantum non-Markovianity. 



II. QUANTUM MARKOV CHAIN 



To define Markov chains, we apply the toolbox containing CP-maps, Kraus matrices, POVMs, selective and non- 
selective quantum measurements, detailed by the monographs 2-6 . We understand by a non-selective Markov chain 
a series of quantum states Qo, Qi, ■ ■ ■ , Qt ■ ■ ■ of a given system where each state along the chain depends but on the 
preceeding state respectively: 

Qt+l=M t Qu t = 0,1,2,..., (1) 

where M. t are CP maps. If not stated otherwise, we restrict ourselves to stationary chains M. t — A CP map can 
always be represented by certain Kraus matrices M\ such that 
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FIG. 2. Selective Markov chain: the refinement of the non-selective one in Fig. 1 1(a)] The post-collisional state of the molecule 
is measured, and yields the random outcome A. 



Accordingly, non-selective Markov chains ([I]) can be decomposed into selective Markov chains given by the recursion 
relation: 

Qt+i = —M xQt Ml (3) 
P\,t 

A may represent the random outcome of a POVM measurement characterized by the matrices M\. Here TvM\Q t M\ 
is the outcome probability. As a matter of fact, the selective chain g t depends on all measurement outcomes A prior 
to t, this dependence is suppressed in our notation. The stochastic average of the selective Markov chain ^ over the 
measurement outcomes A yields the non-selective Markov chain ([I]). 

We can always construct a microscopic model for a given quantum Markov chain. Consider an abstract ideal gas 
(reservoir) of identical molecules each in state £. Independent unitary collisions will generate the CP map M of Q: 

Qt+ i=TT xe8 [U{ et ®t)U\ (4) 



where U is the collision matrix, Tr rcs is the partial trace over the molecule state. The process is shown in Fig. 1(a) 
Any map ^ can be generated by a suitable unitary mechanism Q whereas the choice of U and £ is never unique 
[5J[5]. If we inspect, i.e. measure, the post-collisional state of the molecule we get the underlying microscopic model 
of the selective Markov chain ^ , shown in Fig. [2] As to determining Kraus matrices M\ , their possible choice is 
simple if the molecules are prepared in pure state £ = |\&) 

M A = <A|Z7j*). (5) 

We show concrete examples in Sec. [IVJ 

III. QUANTUM NON-MARKOV CHAIN 



Consider the chain 



Qt=M(t)g ; i = 0,1,2,..., (6) 



where M.(t) is a t-dependent CP map. The chain is Markovian, see (fl]), if A4(t) is divisible, i.e. can be written in the 
following form: 

M(t) = M t M(t - 1); t = 1,2,..., (7) 
with some sequence of CP maps Mi, M2, ■ ■ ■ , A4t, ■ ■ ■ ■ To be clear, we require that 

M{l) = Mx, 

M{2) = M 2 Mi, (8) 
M(3) = M3M2M1, e.t.c. 

If such a sequence Aii, M2, ■ ■ ■ ,-Mt, ■ ■ ■ does not exist, i.e., when A4(t) is not divisible into the same set of factors 
for all t, the chain g t is NM. 

The microscopic mechanism of quantum non-Markovianity can be quite complicated. A recent study [18 considers 
single qubit maps which are not divisible at all. To model such maps microscopically a class of finite NM chains has 
been constructed. Its key mechanism has been independent collisions with entangled molecules, the general structure 



is shown in Fig. 1(c) Another work [19] considered a single qubit environment, interacting unlimited times with 



the system qubit via a controlled-rotation. Our approach will be different and elementary. We construct a class of 
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FIG. 3. Selective NM chain: the refinement of the non-selective one in Fig 
collision) state of the molecule is measured, and yields the random outcome A 



1(b) The post-collisional (i.e.: after the double 
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FIG. 4. Selective one-qubit Markov chain: collisional model. The system qubit g interacts (collides) with the reservoir qubits 
(molecules) £ via XOR-gates; the post-interaction states of the molecular qubits are measured, yielding the random outcomes 
A. 



stationary non-Mar kovianity. Our molecules are independent, and do not become entangled prior to their collisions 
with the system. However, they will collide multiple times — twice, for simplicity's sake — and this constitutes the 



memory mechanism. Fig. 1(b) shows the general NM structure. Unlike the previous NM structure shown in Fig. 1(d) 
ours allows us, by its construction, to insert measurements without changing the physics of the process 
we can insert measurements on each molecule after its second collision, as shown in Fig. [3] 



Obviously, 

Rather than discussing this NM model generally, in what follows we concentrate on its simplest qubit realizations. 



IV. ONE-QUBIT NON-MARKOV CHAINS 



In order to highlight the similarities and differences, first, in Sec. |IV A} we are going to construct a one-qubit Markov 
chain, and then, in Sees. IV B and IV C by extending it, we study two different NM structures. 

Throughout this section the molecules constituting the reservoir are single qubits in pure initial states £ = \^!^} (^<p\, 
where 



= cos0|O) + sin 011) . 
For later convenience, we write the state of the molecule as 

= exp{i<P<r y ) |0) . 



(9) 
(10) 



A. One-qubit Markov chain 



Let our central system be a single qubit, and let us construct a Markov chain ([l][3]), also see Figs. 1(a) and [2] 



Couple the molecule in state £ to the system qubit in state g t via the XOR-gate. The corresponding two-qubit 
unitary operator U is the following: 

U = XOR sys _ mol . (11) 

The XOR-gate installs the unsharp measurement of the system qubit by the molecular qubit. The strength of the 
measurement is controlled by the parameter cj) of the molecular pre-collisional state. If 4> — 0, tt/2, the measurement is 
projective while <j) ss 7r/4 yields a weak measurement |20j : the value 4> = 7r/4 decouples the system from the molecules. 



The two Kraus matrices are defined by M\ = (A| U l^f/^) for A = 0, 1, cf. Eq. (|5|). Inserting (10), we have 



M x = {\\Uexp(i<f>a y )\0) 



(12) 
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To perform the partial trace operation ([i]) from the four-dimensional state U(gt <8 £)£7^ to the two-dimensional state 
g t , we use the deeper level selective form ^ for technical reasons. 

The unitary evolution of the molecule-system compound in the product basis |mol, sys) = {|00) , 1 01) , 1 10) , |11)} 
is given by 



U exp(i0tjj,) 



/cos (j> — sin <f> 

sin cc 
sin <p cos (j) 

\ cos 4> — sin < 



(13) 



The two Kraus matrices in this representation are simply the left block matrices of U: 



V sin <p J V 

The non-selective Markov chain ([I]) is as trivial as [? ] 

e&J 1 = eSo ^i 1 = eii efr 1 = skp^e&i (15) 

i.e., the diagonal elements are preserved while the off-diagonals will step towards zero unless we took the singular 
molecular states with <fi = it/ 4: 



_ (&> o 
g° n 



(16) 



We can say that our Markov chain with 4> =^ 7r/4 is asymptotically equivalent with a single von Neumann projective 
measurement. Even if the single collisions are in the weak measurement regime, their cumulative effect is the projective 
measurement, as is well known, e.g., from [91411] . 



B. Repeated XORs with time overlap 



To construct the simplest non-Markov chain, let each molecule interact with the system qubit twice, and let there 

We are interested in stationary chains hence the 



be an overlap between subsequent collisional periods, see Fig. 5(a) 



pattern of double collision will identically repeat itself along the chain. Note, however, that the starting pattern must 
always be a "broken" one, i.e., in the beginning there is a molecule that collides only once with the central system. 

From the classical theory of NM chains with finite memory time we know that adding suitable amount of memory 
(bits) to the system yields an equivalent Markovian chain. Similarly, we attach a memory qubit to the system qubit. 
We can make e.g., one distinguished molecule the memory. Consider the resulting scheme in Fig. |5(b)| equivalent 



with Fig. 5(a) The system+memory compound is a two-qubit composite system whose state will be denoted by g t 
and, as we see, it undergoes independent collisions with the rest of the molecules of state £ each. Our trick is that we 
repeatedly swap all system-environment entanglement into the system-memory compound. 

Accordingly, we have obtained a Markovian quantum chain for the system+memory compound g t instead of the 
system g t alone. Markovian structure Q-Q applies invariably. The molecule-memory-system compound is a three- 
qubit system. The three-qubit unitary operation U is the following: 



U = XOR, 



sys— mol 



SWAP 



mol— mem 



XOR. 



sys— mol- 



(17) 



In the basis |mol, mem, sys) = {|000) , |001) , 1 010) , . . . |111)}, we calculate the collision matrix (17) times exp(z</>cr y ), 
with notations c = cos 6 and s = sin <b: 



U exp(i(f)tTy) 
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(18) 
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FIG. 5. Equivalent one-qubit NM and two-qubit Markov chains, (a) NM chain: the system qubit g interacts (collides) twice 
with the same reservoir qubits (molecules) £ via two XOR-gates, with overlap between 'nearest' molecules; the post-interaction 
states of the molecular qubits are measured, yielding the random outcomes A. (b) Equivalent Markov chain: the system-satellite 



two-qubit g compound interacts (collides) with the reservoir qubits (molecules) £ via the swap- and two XOR-gates but without 
the overlap between collisional intervals of different molecules; the post-interaction states of the molecular qubits are measured, 
yielding the random outcomes A. 



According to (12), the two Kraus matrices Mo and Mi in the basis |mem, sys) = {|00) , 1 01) , 1 10) , |11)} are given by 
the upper-left and lower-left 4x4 blocks, respectively: 



M 




Mi 




(19) 



Using these Kraus matrices, the non-selective Markov chain ([TJ reads: 



Qt+i 



/c 2 (£oo 

cs (i?30 

cs(qoo 



§22) cs ($>3 

Q12) s2 (e{i 

Q22) S 2 ($)3 

Q12) cs(Qii 



Q21 
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Q21) S Z (Qoo 

£33) cs ($so 



cs ($>0 



022) C 2 (^ 3 

Q12) cs{Qu 

§22) cs ($)3 

Q12) c 2 (q\i 




(20) 



On the r.h.s. we can identify the diagonal part of the system density matrix g 0Q = (qqq + Q22) and Qn = (goo + £33) 
as well as the specific correlation C*_ = Ti(a x ® <J-Qt) = (f?03 + £21) between the memory and the system. Observe 
that all these quantities are invariants along our NM chain. Not surprisingly then, the chain will immediately reach 
the stationary state gi = §2 • • • — 600 fully parametrized by the above invariants: 



csC x+ 

CSQqo 



cs 



X- 



CSQ11 



CSQqq 

s C x+ 
csC x+ 



c 2 C° x - 

CSQu 

csC°_ 
c 2 9°n 



(21) 



The free parameter C x ± depends on the starting "broken" pattern. In the case shown on Figs. 5(a) and 5(b) it 
can be set to zero, assuming uncorrelated initial satellite memory with (a x ) — 0. The above stationary state ol the 
memory+system compound then becomes a separable, disentangled state: 



0oo l**> ® 10) (0| + Q°u \%) (%\ ® |1) (1| 



(22) 
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where the memory state |^) = s |0) + c |1) is not orthogonal to \^^>) = c|0) + s |1) (unless = ir/4). The stationary 
density matrix g± = £ 2 ■ ■ ■ = goo of our central system is trivial: 

Qoo = Tr mom goo = ^0 ^ • (23) 

It is the mixture of the basis states with the initial population, as if go were ideally measured by the collective of the 
molecules, like in case of the Markovian chain, cf. Sec. |IV A| 

The repeated-XOR model has trivial invariants, this is the reason it reaches the asymptotic state in one step already. 



C. Distributed XORs with time overlap 



In order to make a NM chain evolve asymptotically towards the stationary state we must assure that there are 
fewer trivial invariants. Let us go back for a moment to the Markovian chain with single XOR collisions (131 and 
introduce a NM variant. Suppose that each XOR-collision takes a finite time and subsequent XOR operations overlap 
in time. To construct a simplest discrete model of time-distributed XOR, we make XOR in two steps separated by 
unit time. Let both steps correspond to the following square-root of XOR: 



%/XOR 



(24) 



with the 



(Here the convention y/i/2 = e l,r ' 8 /v2 has been chosen.) Our non-Markov model is shown in Fig. 5(a" 
XOR collisions replaced by the above VXOR. 

Obviously, we have the equivalent Markov chain, like in Fig. |5(b)| and we can perform the same calculations as 
before. We calculate the unitary matrix 



U = VXOR, 



■sys— mol 



SWAP mol 

— mem 

VXOR sys 

— mol 



(25) 



times exp(i0(Ty), and read out the two Kraus matrices: 



Mn = 




Mi = 




(26) 



with (3 = e'^/2. These Kraus matrices in ^ yield, after direct calculations, the following Markovian chain |lj for the 
memory-system compound: 



Qt+i- 



(c 2 {qIq + Q22) c/3A t C3($ fl +£>2 2 ) 

cpK t i/2(#i+e|s) s pK t 

cs(goo+i22) s(3A t s 2 ($) +/ 22 ) 

V -c0i&t -2if3 2 (~g[ 1 + Ql 3 ) - S f3iA t 



where 



_ icpA t 
20*(g[ 1 + g 33 ) 

ispAt 
l/2(& + ^3)/ 



(27) 



A t - + gl 3 ) + ($ 3 + £ 21 ) 



and it satisfies a closed evolution equation: 



A t+ i = sm(2<£)A t 



(28) 



(29) 



therefore this is a convenient parametrization. It follows from (27 1 that the matrix elements not containing A will 
take their final stationary values immediately after the first step, just like in our previous model in Sec. IV B However, 
the elements with A show an exponential relaxation (29) toward zero if sin(20) ^ 0. The relaxation of At governs 
the asymptotic diagonalization of the system density matrix: 



gt = Tr„ 



£00° 



(l+«sin(20))A t /2 



(1 - isin(20))A t /2 



t=o 
611 



(30) 



Indeed, for t — oo one has Aoo = 0, the system density matrix g^ becomes diagonal with the initial pop ulations 



£00 1 Qiv As ^° the stationary density matrix g^ of the memory-system compound, the Markov chain (27) yields a 
separable state again, like in our previous model in Sec. |IVB| This time we get 



= g° 00 |**> ® |0) (0| + g° n \%) (%\ ® |1) (1| 



(31) 



|l))/\/2, different from j^) in the model in Sec. 



IV B 



where the memory state |^) = (|0) + i 
to l^^). This is what makes our system quantum NM, a distinction discussed in the next section. 



is never orthogonal 



V. CLASSICAL/QUANTUM NON-MARKOVIANITY 

Our construction of NM quantum chains may become universal if we allow for more than two repeated collisions 
with a single molecule and/or overlaps of collision periods between more than two molecules. Whether or not an 
arbitrary NM chain will be reducible to ours remains an open theoretical issue. Our class of NM chains is unique for 
at least one thing: we can always identify the satellite memory to make the time evolution of the resulting system- 
memory compound a Markovian chain. We can always read-out from the circuit of the NM chain how many qubits 
are needed for the satellite memory. This number of qubits is a useful count of non-Markovianity of our chain, e.g., in 
computational simulation this number gives an upper bound on how much data should be stored dynamically together 
with the data of our system of interest. 

However, this count may be significantly larger than the amount of information to be contained in the satellite. 
Suppose that we have determined the minimum number of qubits needed for the satellite, so the count of non- 
Markovianity is known. Then a suitable informatic quantity, like the mutual information [6] of the satellite memory 
(M) and the system (S) might play the role of non-Markovianity measure /j.nm- I n t ne stationary regime one gets 

mM = I(S : M) = H{S) + H(M) - H(S, M) (32) 

where H stands in turn for the von Neumann entropy of the S-state TrmemPoo = g^, of the M-state Tr S y S £ioo, and of 
the composite state g^. 

We should stop for a second, and distinguish quantum from classical non-Markovianity. The attentive reader may 
have noticed that our construction of quantum NM chains in itself has nothing particular for quantum chains, it roots 
in a similar construction of classical NM chains. In fact, any classical NM chain can be represented by an equivalent 
Markovian chain if we assign a sufficient satellite memory. The minimum size of the satellite memory (e.g.: in bits) is 
the natural count of non-Markovianity, the mutual information may be the measure of non-Markovianity. In the case 
of a generic quantum NM chain, it may happen that the minimum satellite memory still consists of bits, qubits are 
not required at all. In this case we say that the chain is classically NM, and its quantum non-Markovianity is zero. 
It is re marka ble that both NM chains in Sec. |IV| have turned out to be quantum NM. If we look at the composite 
states ( 22|3I I of the system+satellite compound, we see that classical satellite bits would not work, we need satellite 



qubits. Although we got zero stationary entanglement between the system and satellite, it does not mean the lack of 
quantumness. The price of getting rid of quantum non-Markovianity cannot be paid in classical bit instead of qubit. 
The count of quantum non-Markovianity is I for both NM chains. In order to distinguish quantum from classical 
non-Markovianity we use the notion of discord |2I) . The classical measure of non-Markovianity may be defined by 

/4m = J(S : M) {nf } = H(S) - ff(S|{nf }) (33) 

where the rightmost term means the average von Neumann entropy of S when the satellite memory is undergoing the 
projective measurement via the set 11^. The quantum non-Markovianity measure is the discord itself: 

Matm = US : M) - J(S : M) {nf } . (34) 

The sum of quantum and classical non-Markovianity measures yields the total measure Unm- 



We can check our proposal on the two models of Sec. 



IV 



we get non-zero discord. If they were orthogonal, we coulc 



Since the memory states 



'iK ) are non-orthogonal, 



get zero discord, and also we could measure the memory 
after each collision so that a single classical bit 0, 1 could be retained instead of the memory qubit: the chain would 
be classically NM in the stationary regime, with non-Markovianity count I (bit). 




FIG. 6. NM quantum chain of advanced overlap: the system interacts with each independent molecule twice, with overlap 
between collisional intervals of 'nearest' and 'next-to-nearest' molecules. (Boxes connected with vertical lines represent the 
bipartite unitary collisions.) 



VI. BRIEF SUMMARY AND OUTLOOK 



In the framework of the abstract collision model of system-reservoir interactions, we constructed quantum chains 
meeting the minimal requirements for NM stationary time evolution. A method of systematic construction of the 
equivalent Markovian quantum chains using explicit memory allocation is discussed. Due to the transparent NM 
structure, we can always identify a well-defined part of the reservoir as the memory, this part is called the satellite 
memory of the system. The time evolution of the system-plus-satellite-memory compound is Markovian. We suggest 
a novel distinction: a given quantum NM chain is either quantum NM (if the satellite requires qubits) or classical 
NM (if bits suffice for the satellite). Accordingly, we suggest the numbers of satellite qubits and bits, respectively, 
as counts of quantum/classical non-Mar kovianity of a given quantum chain. The mutual information and discord is 
proposed to measure non-Markovianity and quantum non-Markovianity, respectively. 

Two examples of one-qubit NM chains are discussed. In both examples, the corresponding satellite memory is a 
single qubit. The stationary state is exactly calculable, the non-vanishing discord of the satellite qubit w.r.t. the 
system qubit indicates that our examples are quantum non-Markovian. 

Although our calculations are restricted for non-selective NM chains, the structural transparency allows for the 
selective refinements. Option of monitoring the NM system is inherent in the model. 



For achieving a higher level non-Markovianity than discussed in our work (Fig. 1 (b) ) , one can construct a similar 
chain with more collisions per molecule and/or longer intervals of overlap. In Fig. |6j e.g., the number of collisions 
is invariably two whereas the lengths of overlapping intervals have been increased for three collisions. In return, the 
equivalent Markovian structure needs two satellite qubits instead of one, the non-Markovianity count is 2. The ability 
to handle complex NM processes may pave the way to the construction of environmental interaction models which 
better approximate reality. 

Support by the Hungarian Scientific Research Fund under Grant No. 75129 and support by the EU COST Action 
MP100 are acknowledged. The authors thank Tamas Geszti for useful discussions. 
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